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OO ' We consider the scattering of a low-frequency gravitational wave by a massive compact body in 

■ vacuum. We apply partial-wave methods to compute amplitudes for the helicity-conserving and 
helicity-reversing contributions to the cross section, accurate to first order in Muj. Contrary to 
previous claims, we find that the partial-wave cross section agrees with the cross section derived via 

'"^ I perturbation-theory methods. 

in ; 

I. INTRODUCTION 

cr: 

i Gravitational waves - propagating ripples in spacetime - are a key prediction of General Relativity (GR). 

' Despite strong indirect evidence for their existence [l|, Q , they have yet to be measured directly. This is hardly 
surprising given the expected amplitude of waves reaching Earth (with a dimensionless strain oi h ^ 10^'^^). 
, However, nine decades after the formulation of Einstein's theory, many experimentalists are now optimistic 
PsJ ' that "first light" detections are imminent, at either (existing) ground-based [if or (future) space-based 0] 
' interferometers. 

CN ■ 

. In this note, we ask a simple question: when a long-wavelength gravitational wave impinges upon a massive 

compact body, what is the differential scattering cross section? We will assume that the incident wave is 
monochromatic, long-lasting, and sufficiently weak that the gravitational field equations may be linearised. 

■ Thus, the problem is characterised by a single dimensionless parameter, 

> ; Muj^ nrs/X, (1) 

. (with units G — c — 1) which conveniently expresses the ratio of incident wavelength A to the Schwarzschild 

^ , horizon of the compact body. In this paper, we concern ourselves only with the long-wavelength regime, 
■ ■ ■ ' in which Mu ^ 1. 

It is no surprise to find that this simple question has been asked, and answered, by many authors. The 
literature on the gravitational scattering of massless waves of various spin {s — 0, ^, 1, and 2) is extensive 
and stretches back over forty years (see 01 and for summaries). Nonetheless, to our knowledge only one 
paper 0, written in the late 1970s, tackles this problem via partial- wave methods. The purpose of this note 
is to revisit and improve Matzner and Ryan'spioneering study 

Over the years, various authors 0,[i,&E3, U, 12| have shown that, in the long wavelength hmit {Muj <C 1), 



o 



the cross section depends on the spin of the scattered field, as follows: 

s = 0, Scalar wave 
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1 

sin''(e/2) 

cos''(e/2) _ 1 

sin-i(6l/2) — 2 ' 

cos*(e/2) 
sin4(e/2) 

cos (e/2)-|-sin''(e/2) 
81111(6/2) 



Neutrino [bl 

(2) 

1, Photon [c] 

2, Grav. wave [d]. 
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It is worth noting that the gravitational result is somewhat anomalous, in that it doesn't follow the same 
general rule [da/dn = cos^' (6 / 2) / s\n^ {9 / 2)] as the other fields. 

Equation ([2H) is the conclusion of (at least) three separate studies. The first derivation was carried out 
by Westervelt who applied perturbation theory to the linearised gravitational equations. Peters found 
the same result via a Green's function approach, whilst De Logi and Kovacs applied Feynman diagram 
techniques. On the other hand, Matzner and Ryan 0] derived a different and more complicated formula by 
applying partial wave methods. Matzner and Ryan recognised that the lack of agreement between their result 
and the rest of the literature was surprising, given that, in the case of Coulomb scattering, the lowest-order 
partial wave cross section and the Born approximation are in exact agreement. 

The primary aim of this work is to show conclusively that Eq. ([^Ji) is indeed correct by improving the 
partial wave analysis of Matzner and Ryan. A secondary aim is to clarify the origin of the 'extra' term 
sin^(6'/2) in the gravitational cross section ^^). This term is a direct consequence of the non-conservation of 
helicity in gravitational-wave scattering. As we shall see, helicity is not conserved because 'axial' and 'polar' 
waves are scattered ni_ different ways. 

As is well-known 



13j , first-order perturbations to the Schwarzschild metric may be divided into two classes. 



according to their behaviour under spatial inversion. Axial (or odd) perturbations pick up a factor of (—1)'+^ 
under inversion, whereas polar (or even) modes pick up a factor (—1)'. In our discussion, axial/odd modes 
are referred to as having negative parity (P = — 1), whereas polar/even modes are said to have positive parity 
[P — +1). As we shall see, partial waves of the same I but opposite parity pick up different scattering phase 
shifts. 



II. ANALYSIS 



The remainder of this note is organised as follows. In section Hi Al we briefly recap the results of Matzner 
et al. to write the cross section in terms of amplitudes which are expressed as partial wave series. In IIIBl 



we define the phase shifts and discuss their asymptotic values IJ] in the long-wavelength regime [Mlo <^ 1 



In III C[ we derive a useful formula for the spin- weighted spherical harmonics. In IIID( we show that, in the 
long-wavelength limit, the partial wave series may be summed to give pleasingly simple results. We conclude 
with a brief discussion in section IIIII 



A. Partial Wave Series 

In the late 1970s, Matzner and co-workers 0, [l5, l3] showed that the differential cross section for the 
scattering of gravitational waves from a spherically-symmetric compact object can be written as the sum of 
the square magnitude of two amplitudes, 

^ = |/(x)p-f|g(x)p. (3) 
These amplitudes may be expressed as partial-wave series, 

oo 

fix) --EE bM-^^Sf,) 1] -2Yi'{l) -2Y,\x), (4) 

OO 

In these expressions, exp(2i(5;^) are phase factors to be determined from a radial equation, sY[^{x) are spin- 
weighted spherical harmonics, and x = cos 0, where 9 is the scattering angle. Note the presence of the sum 
over even and odd parities, P — ±1. 
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The first amplitude f{x) corresponds to (that part of) the interaction which preserves the helicity (i.e. 
for wliicli the helicity of the scattered wave is the same as the helicity of the incident wave). The second 
amplitude g{x) corresponds to (that part of) the interaction which reverses the incident helicity. As we see in 
the next section, the helicity-reversing amplitude is non-zero, because the phase shifts Sf^ depend on parity 
P. In this respect, gravitational wave scattering is unlike scalar, neutrino, or electromagnetic scattering. 



B. Phase Shifts 

The phase shifts of odd parity (P = — 1) may be found from Regge and Wheeler's radial equation, 

+ 6M 



which describes axial perturbations. Here, is a tortoise coordinate defined by dr/dr^ = 1 — 2M/r. To 
find the phase shifts, one must solve this equation subject to the ingoing boundary condition at the horizon 
[R{r) ~ exp(— iwr*) as r — > 2M]. The asymptotic solution in the far-field is 

i?(r) ^ Ain e-'^"-- + Aout e*""-* , as r ^ oo, (7) 

and the phase shifts are determined by the ratio of the ingoing and outgoing coefficients, 

eM'^^SfJ = (-l)'+lAoutMin. (8) 

The phase shifts of even parity (P = +1) may be found by solving Zerilli's radial equation [3]. It has been 



shown [5] that the phase shifts of even parity are related to those of odd parity by 

exp(2z5+ ) - i^ + m + ml~l) + l2^Mc. 

Thirty years ago, Matzner and Ryan Q conducted a partial wave analysis in the low- frequency limit. They 
assumed the odd-parity phase shift to be approximately 



lim^exp(2z<5,„)«^,)^^^^-^. (10) 



r(/ 1 - 2iMuj) 
a/^^'^^o""^'"""'"' r(Z + 1 + 2iMuj) 

More recently, Poisson and Sasaki showed that the exact result for the phase shift in this regime is actually 



where 



lim^exp(2i(5,";J = e"** exp (-4iAfw/3/) , (11) 



A = U*(^ + i) + *(0 + ^^7rr^)' ^ith *(0^iin(r(0). (12) 



2 V /(; + i) 

Here, $ = — 4Afti> ln(4Ma;) is an overall phase factor which has no effect upon the cross section. The result 
of Poisson and Sasaki can be rewritten 

exp(2i5;;„) = e-'*-2»^^" exp [^AiMuj'fil + 1) + 8iMuj/l{l + 1)] + 0{M^uj^) (13) 

„ -i'S>-2^Mi^ SiM<^/lil + l) r(^ + 1 - 2lMuj) »,,-2 2^ /, .^ 

r{i + i + 2iMu) ^ ^ ^ ^ ^ 

The extra factor of pj-gg^j^^ p4p but not in PH)) proves significant, as we see in section Hi Dl 

But first, let us briefly digress to study the spin-weighted spherical harmonics. 
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C. Spin- Weighted Spherical Harmonics 

To compute the amplitudes ([4]) and ([5]) we require expressions for the spin-weighted spherical harmonics 



-2Yi^{x). These may be found by acting on spherical harmonics of spin-weight zero, oYi^{x) = y ^^^i Pi[x), 
with ladder operators [l^. The spin- weight is lowered with the operator S, and the azimuthal number is 
raised with L+. These operators are defined by 

SsYrix) = (^v/T^9. - ^i^) sYrix) = -vii + sKi-s+i)s-iYr{x), (is) 

L+,Yrix) - - (^Vl-x^d, + sYrix) = ^il-m)il + m+l),Yr+\x). (16) 

Here, dx is shorthand for the partial derivative with respect to a; = cos 6'. By acting with SL^SL^ on qY^^Ix), 
it is straightforward to show that the spin-weighted harmonics in ^ and ^ can be written 



-2^/ (^) - V il-l)l{l + l)il + 2) ) ■ 

Their values in the forward and backward directions are particularly simple. 



D. Scattering Amplitudes 

In this section we show that, to first order in A'luj, the scattering amplitudes f{x) and g{x) are given by 

^ ^ ^(l + 2^M^) [sin2((?/2)] 

Af-i g(6») = e"^* X sin2(6'/2) + 0(Mw), (20) 

where e~** is an irrelevant phase factor. 

First, let us consider the helicity-conserving amplitude / defined by Using results ([H]) . and PT|) . 
it may be written as 



/(:e) = e-'*-^^^'^-(l + x)^ a,(l - x)dAil - x)dxF{x), (21) 



where 



rjx)- ' V (2^ + l)e^'^'P.(x) 

^ ^" 2*c^A^ (/-l)U; + l)(^ + 2)' ^ ^ 

and 



/=2 



- [{1-1)1(1 + + 2)- 12iMu;J T{l + l + 2iMLo) +U[iuuj). [Z6) 

To find the amplitudes in the long- wavelength limit, it is only necessary to keep terms up to first order in 
AIlu. Hence 

il-mi + l){l + 2) _ (/__l)Z(; + l)(Z + 2) ^^^^ 



(/ - 1)1(1 + l)(l + 2) - UiMuj (l-l- 2iMoj)(l + QiMuj)(l -I- 1 - QiMoj)(l + 2 + 2?;Ma;)' 

p8.Mc./;o+i) ^ (l + ^iMuj)(l + l~%iMuj) 
and e . (25) 
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The phase factor e ' can then be written 

The series F{x) may now be computed with the aid of Eq. (7.127) from Gradshteyn & Ryzhik 20|, which 
gives 

To find F{x) we substitute a = 1 + 2iMuj into the above formula and compare with J_^F{x)Pi{x)dx. This 
imphes that 



2-(l + 2»*=f") T(l-2iMuj), , . 

Fix) = M- U\ - 28 

^ ' (2iMw)2(l + 2iMa;)2r(l + 2iMw)^ ' ^ ' 

Plugging ([28|) into (f2T|) and taking four derivatives yields the pleasingly simple result 



,r(l-2iMcj) ■ 
r(l + 2iMa;) [i (i _ ^.y ^-^^^"^ 



/(X) = Me-— lii^^ ..,^^":a.M. (29) 



which is the same as (|T9|). 

The helicity-reversing amplitude 5(2;) defined in ([5]) may be written 

. £:;^(12.M., ^ _,r=(-., (30, 

The denominator of this expression ensures that the series converges quickly. We are only interested in the 
amplitude to first order in Afw, and there is already a factor of YliMuj in ([50]) . Hence it is justified to take 
g-«Af<^A ~ 1 + C'(Mw). 

In order to compute g{x)^ let us first consider the integral 



^ ^ ' {\-xfd^{\ + x)dM^ + x)d^Pi{x). (31) 



4^ (/-l)Z(/ + l)(? + 2)7_i 
Integrating by parts four times we obtain 



-OA + l (-1)' 

"^*V 47r (/-l)Z(Z + l)(Z + 2) ^^^^ 
for I > 2. Applying the orthogonality relation for the spin- weighted spherical harmonics, 
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-2Yi'ix)_2Yif{x)dx^—Su', (33) 

I ZTT 

we conclude that 

L '^^^ (-^^'^ ^ ''''^'''\IW ii-ml{)ii + 2y (^4) 

Hence, by comparison with ([5^ . 

ff(x) = iMe-^*(l-x), (35) 

which is the same as result 
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III. DISCUSSION AND CONCLUSION 

In the preceding sections we have computed the partial-wave scattering amphtudes defined by Matzncr and 
co-workers 0, EB] in the long- wavelength limit {Mlo <^ 1). Substituting p9| and (pO)) into ([3]) we conclude 



that 

dcr _ 2 sin* (61/2) -t- cos** (61/2) 
sin4(0/2) 

The partial- wave result is therefore consistent with all previous studies \\\ , which reach the same result 
by applying perturbation-theory methods. 

It is worth remarking that, unlike other polarised waves (i.e. neutrino and photon waves), the scattering 
cross section of the gravitational wave ([36|) is non-zero in the backward direction (Q — -k). As we have seen, this 
is because the phase shifts in the partial wave series are parity-dependent (Eq. [9]) . This implies the existence 
of a helicity-reversing amplitude gix). As noted by De Logi and Kovacs [11], "if the incident radiation is in 
a pure helicity state, the backscattered {6 = tt) radiation must have the opposite helicity". 

Finally, we note that a spherically-symmetric interaction will not induce a net polarisation in an initially 
unpolarised beam. Right- and left-circular polarisations are scattered in the same way (Eq. I36|) . If the 



scattering body is rotating, the spherical symmetry is broken. It has been suggested ll| that, to first order 
in Mu, polarisation is not induced by a rotating scatterer. However, outside this regime (i.e. when Muj ^ 1) 
we would definitely expect to see some polarisation effects caused by rotation. For instance, if the incident 
wave were to impinge along the rotation axis of a Kerr black hole, then the co-rotating polarisation would 
be enhanced by the superradiance effect [l3| • We hope to make a numerical invesigation of scattering from a 
rotating black hole in the near future. 
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